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Abstract 

We present a general numerical method for computing precisely the false vacuum decay rate, 
including the prefactor due to quantum fluctuations about the classical bounce solution, in a self- 
interacting scalar field theory modeling the process of nucleation in four dimensional spacetime. 
This technique does not rely on the thin-wall approximation. The method is based on the Gelfand- 
Yaglom approach to determinants of differential operators, suitably extended to higher dimensions 
using angular momentum cutoff regularization. A related approach has been discussed recently by 
Baacke and Lavrelashvili, but we implement the regularization and renormalization in a different 
manner, and compare directly with analytic computations made in the thin- wall approximation. 
We also derive a simple new formula for the zero mode contribution to the fluctuation prefactor, 
expressed entirely in terms of the asymptotic behavior of the classical bounce solution. 
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I. INTRODUCTION 



The phenomenon of nucleation drives first order phase transitions in many applications 
in physics, most notably in particle physics, condensed matter physics, quantum field the- 
ory and cosmology. The semiclassical analysis of the rate of such a nucleation process was 
pioneered by Langer who identified a semiclassical saddle point solution that gives the 
dominant exponential contribution to the rate, with a prefactor to the exponential given 
by the quantum fluctuations about this classical solution. The nucleation rate is given by 
the quantum mechanical rate of decay of a metastable "false" vacuum, </>_, into the "true" 
vacuum, </>+. Decay proceeds bv the nucleation of expanding bubbles of true vacuum within 
the metastable false vacuum This picture may also be extended 

to finite temperature field theory [3. llC |. Computing the semiclassical prefactor requires 
the computation of the determinant of the differential operator associated with quantum 
fluctuations about the classical solution. This is a technically difficult problem. Thus, it 
is not possible to compute the fluctuation prefactor analytically, and so various approxima- 
tion techniques have been developed. The most widely-studied is the so-called "thin- wall" 
approximation, in which the bubble wall thickness is small compared to the bubble radius 
SI Q. This provides an ele g a„ t pn y8 ica, picture of nucleation and in this limit certain parts 
of the decay rate computation may be done analytically. A limited amount is also known 



analytically concerning the expansion beyond the leading thin- wall limit, in two dimensions 

nnn nn n 

p, llH, H2I , three dimensions [13, H4| , and four dimensions (151 • Other important approaches 



tensions jlolj . utner important approg 
the average effective action 



to metastable decay use direct lattice simulations [16j . the average effective action [Tj|, and 
phase shift techniques jis| . 

The point of this present paper is to reduce the calculation of the decay rate to a straight- 
forward numerical computation, without relying on any such approximation or expansion. 
The starting point for our approach is an extremely elegant and simple method, due to 
Gelfand and Yaglom for computing the determinant of a one- dimensional differential 
operator. This technique is ideal for numerical implementation. However, its naive gen- 
eralization to higher dimensions is divergent jj^. This is because in higher dimensions 
renormalization is important, and so we must regulate and renormalize the determinant. 
We present an analytic method of doing this using an angular momentum cutoff, and apply 
it to the problem of false vacuum decay. (This method has previously been used to compute 
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the quark mass dependence of the fermion determinant for quarks in the presence of an 



instanton background 



2l|.) A related approach to false vacuum decay, also based on the 



221, and 



Gelfand-Yaglom formula, has been developed recently by Baacke and Lavrelashvili 
we comment in Section II V Bl on the similarities and differences between our approaches. 
Consider the Euclidean classical action 



d*x[-(dM 2 + U{<f>) 



(1-1) 



where the potential U(<f>) has two nondegenerate classical minima, 0_ < <f) + , with £/(</>_) > 
U For quantitative computations, and comparison with previous work, we consider the 



standard quartic potential [2J, |4, 
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whose form is illustrated in Figured 



;i-2) 



The parameter e represents a constant external source breaking the degeneracy of the double 
well potential. In four dimensional spacetime, the mass dimensions of the couplings are: 

U(0) 




FIG. 1: Field potential U{4>) showing the true and false vacua, <p + and <f>—, respectively. 



[A] = 0, [a] = 1, and [e] = 4. We choose a > 0, A > and e > 0, in which case the two 



minima are <fi± = ±a(l ± ^4 +•••)• Note that 

U(4>„)-U(4> + ) = e\l + 



Xa 4 



1.3) 



so for e Aa 4 , we see that e has the physical interpretation as the potential energy difference 
between the two classical vacua. This small e limit is known as the "thin-wall" limit \M 
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because in this case the bubbles of true vacuum within the false vacuum have thin walls 
compared to their radius. 

Expanding the field about the false vacuum 0_ 

= 0- + ^ , (1-4) 

and keeping terms up to dimension four, we find the potential, which is often considered 
directly in the literature j^] 

f%) = ^V- W 3 + ^ 4 . (1.5) 

Here m 2 and rj are related to the original couplings by 

m 2 = ^(30 2 -a 2 ) ; V = . (1.6) 

In order to describe semiclassical tunneling, it is useful to rescale the field ip and the space- 
time coordinates as 

2 

m 

x = mx ; lo = — $ (1.7) 
Then the classical action in terms of these dimensionless quantities is : 



~(W + ~* 2 -~* 3 + f* 4 



;i.8) 



The overall factor is determined by the dimensionless parameter 



m 2 1 ( e \ . , 

"=4?=a( 1 -2^ + -') ' <- h9) 

For our semiclassical tunneling analysis, we assume ,8 > 1. The quartic coupling strength 
in ()1.8|) is determined by the dimensionless quantity 

6 ft i n\ 

a = l^ = 1 -2^ + - ' (L10) 

which tends to 1 in the "thin- wall" limit where e — > 0. The dimensionless parameter a 
determines the shape of the potential, and its deviation from 1 is a measure of the vacuum 
energy difference relative to the barrier height. Figure 121 shows some plots, for various values 
of a, of the dimensionless potential 

1 1 ru 

W) = 2 $2 "2 $3+ 8 $4 ' (L11) 



U(S) 




The false vacuum decay rate per unit volume and unit time is denoted 7, and its one-loop 

□□□□□□□□ 



expression is 
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2tt 



det'(-n + tT($ci)) 



-1/2 



det (-□+[/"($_)) 

where the prime on the determinant means that the zero modes (corresponding to trans- 
lational invariance) are removed. Here $ c i is a classical solution known as the "bounce" 
solution [4] , defined below, and the prefactor terms in (jl.l2j) correspond to quantum fluctu- 
ations about this bounce solution. The second term in the exponent, t^tS^d], denotes the 
counterterms needed for renormalization. The computational challenge is to evaluate the 
rate 7 given a particular form of the classical potential. For our particular quartic model 
this corresponds to computing the rate for various values of (3 and a, the dimensionless 
parameters defined above in (11.91) and (11.101) . 

In the language of quantum field theory |3, |4|, |2J, 7 is half the imaginary part of the 
generating functional of the connected Green's functions with the constant external source, 
— ^, divided by the spacetime volume factor, and so we are essentially computing the 
renormalized effective action for this system [2J, |24|. It is not possible to compute this 
tunneling rate 7 analytically. Indeed, it is not even possible to find the classical bounce 
solution, $ c i, analytically, let alone the quantum fluctuations about this classical solution. 
We present here a simple numerical technique to compute 7 for general forms of U(<p). This 



technique involves a combination of an analytical computation and a numerical computation. 
The analytical part of the computation is related to the regularization and renormalization 
of the one loop effective action. The numerical part is elementary, and can be implemented 
straightforwardly and efficiently in Mathematica. 

In Section IHl we describe how to compute the classical bounce solution numerically with 
very high precision. In Section IIHI we describe how to compute the determinant prefactor 
arising from quantum fluctuations about the classical bounce solution. In Section IIVI we 
explain how to regularize and renormalize our answer using the angular momentum cutoff 

n 

regularization and renormalization scheme developed previously in |2l| . Our approach is 
related to an elegant technique presented recently by Baacke and Lavrelashvili 22], and in 
Section |TV| we also compare and contrast these two methods. In Section |V| we conclude with 
some general comments about possible further applications, and in the Appendix we present 
the derivation of a simple new formula for the contribution to the decay rate 7 coming from 
the zero modes. 



II. COMPUTING THE CLASSICAL BOUNCE SOLUTION 

The first step in computing the false vacuum decay rate 7 is to find the classical bounce 
solution, $ c i(r), which is an 0(4)-symmetric stationary point of the classical Euclidean action, 
with $ c i(r) interpolating between the false and true vacuum as r goes from to 00 jj, l^ . 
The action evaluated on the classical bounce solution determines the leading exponential 
factor in ()1.12|) . and the quantum fluctuations about the classical bounce solution lead to 
the determinant prefactors in (|1.12|) . 

The bounce $ c i(r) solves the nonlinear ordinary differential equation 

-^,-^1 + ^1-^ + ^1 = (2.1) 

with boundary conditions 

$' cl (0) = (2.2) 

$ c i(r) -> = , as r -> 00 . (2.3) 

It is not known how to find $ c i(^) analytically in any nontrivial field theory. Much work 
has been done in the thin-wall approximation, which in practice means finding $ c i as an 
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expansion about the point a = 1, where the two vacua are degenerate. The point of 
this present paper is to reduce the calculation of the tunneling rate 7 to a straightforward 
numerical computation, without relying on any such approximation or expansion. Thus, we 
begin by determining $ c i(r) numerically. 

To obtain extremely good precision we use a combination of both forward and backward 
shooting to compute $ c i(?")- in forward shooting we numerically integrate (|2.1|) [using 4th 
order Runge-Kutta], starting at r = 0, and we adjust the initial value $0 = $ci(0) until the 
second boundary condition ()2.3|) is satisfied to a certain level of accuracy. Since one cannot 
start exactly at r = 0, we begin at some very small r = ro and use a Taylor expansion to 
write 

Mn>) = $ + ^rl (2$ - 3$g + a$jj) 

*d(n>) = \ro (2$ - 3$jj + a& ) . (2.4) 

3 ci 
4 

3 

2 

1 
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FIG. 3: Plots of the bounce solution <3? c i(r) for various values of a: a = 0.5, 0.9, 0.95, 0.96, 0.97, 
0.98, 0.99, with the plateau ending farther to the right for increasing a. Observe that as a — > 1, 
the sharp falloff in $ c i occurs at r ~ j^- , and $ c i can be approximated by a step function. 

We adjust the parameter $ until the large r boundary condition (j2.3|) is satisfied. In 
backward shooting we begin the numerical integration [also using 4th order Runge-Kutta] 
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at some very large r = R, with starting conditions 

K X {R) 



K 2 (R) 



(2.5) 



and adjust the parameter <3>oo until the r = boundary condition ()2.2j) is satisfied. 

We first use forward shooting, which produces a good estimate of <3>o- Then using this 
bounce solution we estimate $oo, and use this as a starting point for shooting backwards, 
which further refines this value of $oo and also leads to a refined value of $o- It is a simple 
matter to obtain 20 - 30 decimal precision for each of $oo and $o m this way. Some bounce 
profiles for various values of the coupling parameter a are shown in Figure H3 
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FIG. 4: Plots of the fluctuation potential U" ($ c \(r)) for various values of a : a = 0.5, 0.9, 0.95, 
0.96, 0.97, 0.98, 0.99, with the binding well of the potential appearing farther to the right for 
increasing a. Observe that as a — ► 1, the potential C/"($ c i(r)) is localized at r ~ j^— , and is 
approximated well by the analytic form in (|2.7|) . 

Given the bounce solution, $ c iO"), the corresponding radial fluctuation potential is 

U"^ cl (r)) = 1 - 3* d (r) + ^(r) . (2.6) 

Clearly, since $ c i(?") is only known numerically, the fluctuation potential U"(<& c \(r)) is also 
only known numerically. Figure 0] shows some profiles of this fluctuation potential, corre- 
sponding to the various bounce profiles in Figure El Note that this fluctuation potential 
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is highly localized, with the localization radius depending strongly on a. As a — > 1, the 
fluctuation potential binds states at radius r ~ the radius of the bubble of false vac- 
uum. Moreover, in this thin- wall limit, the minimum of the potential approaches —1/2, and 
outside the binding well it approaches 1. In fact, as noted in Langer's original work |l[ 
this limit C/"($ c i(r)) is approximated well by the analytic potential 

1 



in 



1 — - sech 2 ( - 
2 I 2 



1 — a 



(2.7) 



III. COMPUTING THE DETERMINANT PREFACTOR 



Since the bounce solution $ c i(r) is a function of r, the fluctuation operator [—□ + [/"( $ cl )] 
can be decomposed into partial waves, with (dimensionless) radial operators 

d 2 3d 1(1 + 2) 



M 



(0 



+ 1 + V(r) 



(3.1) 



dr 2 r dr r 2 

of degeneracy (I + l) 2 , with I = 0, 1,2, . . .. Here the radial potential V(r) is equal to the 
fluctuation potential (j2.6J) with its large radius asymptotic value, 1, subtracted: 

3a 



\Z(r) = -3$ cl (r) + ^$ c 2 1 (r) 



(3.2) 



Likewise, the free fluctuation operator [— □ + £/"($_)] can be decomposed into partial waves, 
with radial operators 

d 2 3d 1(1 + 2) 



Affrcc 



dr 2 



r dr 



+ 



(3.3) 



also of degeneracy (/ + l) 2 . 

For each /, the ratio of the determinants of M.q\ and - M^ e can be computed efficiently 
and precisely using the Gelfand-Yaglom method 



20, 2t 
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28, 



29 



This result 



states that for radial operators of the form (|3. 113.3(1 . 



det(M(i)) 
det(M^ c ) 



lim 



(l+l) 2 



(3.4) 



Here ip^ and ip^ e are regular solutions of 



M(i) ipQ) 



(3.5) 



with the same leading behavior at r = 0. By regularity we can choose this small r behavior 
to be 



~ r l , r ->• 

v$r ~ r ' > r - • ( 3 - 6 ) 

This normalization choice fixes the free solution to be 

i/>$*(r) = 2 l+1 {l + (3.7) 

where 7j + j(r) is the modified Bessel function. 

In fact, a numerical improvement comes from realizing that both i/}(i)(r) and ^>£? e (r) firow 
exponentially fast at large r, so it is numerically better to integrate directly the ratio 



22 



which satisfies the equation 



, (3.8) 



with simple initial value boundary conditions: 

T (0 (0) = 1 ; T ( ' (0) = . (3.10) 
Then the result (13. 4 J) is recast as 



2 



^H = (T (0 (oo)) ( ' +1) - . (3.11) 
det(.M^ c ) un /; v 7 

We stress that the result (J3.ll)) provides a remarkably simple technique for computing 
the determinant of a radial differential operator. It does not require any detailed knowledge 
of the spectrum of the operator whose determinant is being computed, nor does it require 
evaluating and numerically integrating the associated phase shift. Furthermore, the result 
(13.11)1 is ideally suited to numerical evaluation, as initial value boundary conditions are 
straightforward to implement numerically. 

There are three different types of eigenvalue of the fluctuation operator, each having a 
different role physically and mathematically. 
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Negative Eigenvalue Mode : (1 = 0) 

The lowest eigenvalue mode in the I = sector is a negative eigenvalue mode of the 
fluctuation operator, and is responsible for the instability leading to decay. It can 
be identified with homogeneous swelling and shrinking of the bubble of true vacuum. 
This mode contributes a factor to the decay rate_7 related to the absolute value of the 
determinant of the I = fluctuation operator Q, ^ . This determinant is computed by 
numerically integrating ()3.9|) for I = : 

-1/2 



detM 



(1=0) 



detMj- 0) 



|T (0) (oo)P 1/2 . (3.12) 



Zero Eigenvalue Modes : (/ = 1) 

In the I = 1 sector there is a four-fold degenerate zero eigenvalue of the fluctuation 
operator. Physically, these four zero modes are the Goldstone modes associated with 
the breaking of translational invariance. Integrating over the corresponding collective 



coordinates 



produces the factors of |^ in ()1.12j) . In computing the rate 7, we need 



elq. 



We 



the determinant of the fluctuation operator with the zero modes removed 
have found the following simple new formula for the I = 1 sector prefactor contribution 
(see Appendix A for a derivation): 



gcibd] V / det'M {l=1) 
2tt J \ detMf™, 



" 1/2 r ( 3 xv:> 



(3.13) 



where $0 = ^ci(O) is the bounce solution evaluated at the origin, and $oo is the 
coefficient of K\{r) jr in the large r behavior of & c \(r), as in ()2.5j) . The advantage of the 
result (j3.13|) is that it is expressed entirely in terms of the asymptotic behavior of the 
classical bounce solution $ c i(r). This asymptotic information is already generated in 
the precise numerical determination of the bounce solution, so no further computation 
is needed to extract the zero mode contribution to the prefactor. 

Positive Eigenvalue Modes : (/ > 2) 

For / > 2, the fluctuation operator has positive eigenvalues, each of degeneracy (Z + l) 2 . 
These modes correspond to deformations of the bubble shape and thickness. For each 
/, the associated radial determinant is computed simply by numerical integration of 
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FIG. 5: Plots of Tm(r) for Z = 0, I = 1, and I = 2. These plots are for a = 0.5. Note that the 
asymptotic value, TW\(oo), is negative for / = 0, zero for I = 1, and positive for / = 2, confirming 
the discussion in the text concerning the three different types of modes. 



()3.9j) with the initial value boundary conditions ()3.10j) : 

-1/2 

= [ T (0(°° 



detA* 



(0 



(/+l) 2 /2 



(3.14) 



Figure El shows plots of T^(r) as a function of r, for 1 = 0,1 = 1, and I = 2. These plots 
are for a = 0.5. Observe that T(j)(r) approaches very rapidly its asymptotic value, T^)(oo), 
for r outside the range of the binding well of the fluctuation potential (compare with Figure 
H]for ol = 0.5). Also note that the asymptotic value, T^(oo), is negative for I = 0, zero for 
1 = 1, and positive for 1 = 2, confirming the discussion above of the three different types of 
modes. 

For each partial wave with I > 2, the radial determinant, [T(j)(oo)] , is finite and 
simple to evaluate. In the discussion of renormalization it proves convenient to consider the 
logarithm of the determinant factors appearing in the rate (jl.l2j) : 



1, / detM 
In 



2 " I detA^{^ c 



-i(/ + l) 2 lnT w (oo) 



(3.15) 



For large I we can use the radial WKB analysis of 
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lnT( Z )(oo) ~ 



(1 + 1) 



drr V(r) 



32j to find the leading behavior: 
1 



O 



(3.16) 
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In T (1) (co) 




FIG. 6: The I dependence of lnT(n(oo) for three different values of a: a = 0.01 (top curve), 
a = 0.7 (middle curve), and a = 0.9 (bottom curve). The dots show the values of lnT(n(oo) 
evaluated numerically using (|3.9|) . while the solid lines show the WKB prediction Q3.16JI of the 
leading large I behavior. Notice the excellent agreement for I > 20. 

The leading large / behavior ()3.16j) is illustrated in Figure El for several different values of 
the quartic coupling parameter a. Notice the very good agreement between the numerical 
results for lnTm(oo) [solid points] and the WKB estimate [curves from (|3.16j) ]. Further 
subleading corrections to ()3.16|) are computed below in the next Section - see (|4.14|) . 



IV. REGULARIZATION AND RENORMALIZATION 



The large I behavior ()3.16|) means that the formal sum of contributions ()3.15|) to — In 7, 

00 

-^(/ + l) 2 ln[T (0 (oo)] , (4.1) 



2 

1=2 

is quadratically divergent. This divergence should not be too surprising, as we have neither 
regulated nor renormalized the determinant prefactor in the expression ()1.12j) for 7. In the 
language of quantum field theory, we need to compute the renormalized one-loop effective 
action for this interacting scalar field theory U, 0, EJ. Here we will apply the angu- 
lar momentum cutoff regularization and renormalization technique developed in 2l( which 



was successfully applied to the computation of the quark mass dependence of the fermion 
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determinant in an instanton background in QCD. 



A. Regularization 



The first step is to introduce a regulator mass /x using the proper-time representation of 
the logarithm : 



In 



det M 



det Mf^ c 



ds 



{fi 2 s) £ tr 



(4.2) 



rcg 



where the space-time dimension is extended to 4 — 2e, and we have explicitly extracted the 
degeneracy factor (I + l) 2 from the trace of the radial operators and defined in 

f)3.1)l and (jHISl)- We then split the sum over I of these regulated terms into two parts: 



1=2 



In 



detM 



(i) 



detM f fi c 



rcg 



1 L 

^E 



z=o 



111 



det M 



(i) 



detMg c 



1 oo 

+ 5 E 



rcg 



l=L+l 



hi 



det 



(0 



det M^ c 



= Ei 



(4.3) 



Here L is a large but finite integer (in practice, we take L to be of the order of 50 to 100, 
depending on the value of a). Since Si is a finite sum we can safely remove the regulator /j 
and write 

L 



Ei = - 
2 







T (0 )(oo) - 2 In 


**-( 



a 



i^(/ + l) 2 lnT (0 (oo) (4.4) 



i=2 



The first term in Si corresponds to the negative mode contribution (|3.12|) . the second term 
is from the zero mode contribution (|3.13|) . and the final sum corresponds to the contributions 
in partial waves with 2 < I < L. Each term T^(oo) in this last sum is evaluated numerically 
by solving (|3.9|) . The second sum in ()4.3|) 



1 r°° We 



i=L+l 



(4.5) 



cannot be computed numerically using ()3.9|1 because of the infinite sum and the presence 
of the regulator mass \i. Instead, we use radial WKB, which is a good approximation for 
large I, to compute analytically the leading large L behavior of S 2 . This is a straightforward 
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computation using the results of |2l|, |32| • The result is: 



1 r_ (l+ikl + 2) r drrV(r) + ^L r drr3v(v+2) 

2 I 4 Jo 8 Jo 

~T6 r dr (^ + 2 -7i? + ln^) r 3 \/(\/ + 2)|+oQ] . (4.6) 



o 

Here 7^ is Euler's constant. Note that this WKB computation reveals large L divergences 
going like L 2 , L and InL, in addition to a term which is finite in the large L limit. These 
last two terms are exactly of the form expected for renormalization, as is discussed in the 
next section. 

The most important observation at this stage is that the large L divergences of £2 cancel 
exactly the large L divergence found numerically in £1, leaving a finite answer. Indeed, com- 
paring ()4.6|) with (j3.16|) we see immediately that the quadratic divergence cancels. Extending 
(I3.16J) to include 2nd order WKB contributions proves the cancellation of the sub-leading 
divergences also. 

B. Renormalization 

In standard perturbative renormalization theory, the self-interacting scalar field theory 
described by the action (jl.lj) involves two renormalization parameters, A and a 2 , together 
with the wave function renormalization Z. The renormalization counter terms can be found 
by replacing those parameters with A + 5X, and a 2 + 8a 2 . The one-loop approximations can 
be found in quantum field theory text books (see, e.g. J3|) : 

and Z = 1, in the MS renormalization scheme using dimensional regularization. With this 
replacement we identify the renormalization counterterms as 



5 ct S = I d 4 x 



5X i 2 2\ 2 ^ a 



2 



1 r°° (\ 

1 • 1 1 „3i 



(4.9) 

32 i ^ { 1 e -^ E ) r ' iV ( V + 2 ^ ( 41 °) 
In obtaining (PTtTIJ) from flOJ) we have used ((131), CC3)> dH3)> (HHH) , and to translate 
from the original physical field (ft, and couplings A and a 2 , to the field $ and coupling a 
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which correspond to expanding the field about the false vacuum 0_. Notice the appearance in 
(j4.10j) of the combination V(V + 2) when the counterterm is expressed in terns of these fields. 
Furthermore, note that the counterterm only involves the pole-term - and 7^, independent 
of any other parameters, for instance the mass parameter, in the MS renormalization scheme. 
We identify this counterterm (j4.1())) within the WKB result ()4.6j) . with precisely the right 
coefficient and structure. We choose this particular MS renormalization prescription in order 
to compare with the work of Baacke and Lavrelashvili [22j , who also use an MS prescription. 

Combining T,i and S 2 with the above renormalization counter term 5 ct S, we get the 
one-loop effective action: 



MS 



-In T, 



(0) 



00 



2 In 



2 $00 



3,, a n 



1 1=2 



drr 3 V(V + 2) 



1 

16 



drrV{r) 

dr [ I -i In-) r 3 V(V + 2) (4.11) 



Here, the first two terms correspond to the I = contribution (for the negative mode), 
computed using ()3.12|) . and the I = 1 contribution (for the zero modes), computed using 
f!3.13|) . The renormalized effective action Ty^ is finite, and we obtain excellent convergence 
(accelerated by third order Richardson extrapolation j^) for L of the order of 50 to 100, 
depending on the value of a. This result for is plotted in Figure [3 as a function of a. 

We can alternatively express (j4.11j) with the L dependent WKB subtraction terms in- 
cluded as subtractions inside the I sum : 



MS 



-In T, 



(0) 



,00 



2 In 



— <J> 

2 00 



a 



1=2 



)) 



\ j °° dr rVir^ + \j™ dr r 3 V(V + 2) 



(l + l) (H 

3 f°° I f°° /l r 

■-J drrV(r) + —J^ dr r 3 V(V + 2) (- - lE - In - 



(4.12) 



Here we have simply used 

L 



E( / + 1 ) = ^ L + 1 )( L + 2 )- 3 ; E 



InL 



IE 



(4.13) 



1=2 1=2 
Thus, (14.12)) shows that we can extend to next-to-leading order the large I behavior in (J3.16|) 
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FIG. 7: Plot of the MS renormalized effective action 1)4.11)1 as a function of a. 
of In (T(z)(oo)), the logarithm of the determinant of the radial operators in (13.1)) : 

We have confirmed numerically that the remainder is indeed O (jj^s) ■ Note that 1)4.14)1 
provides a simple closed-form expression for the large I behavior of In (T(y(oo)), and this 
expression is local in the fluctuation potential V(r). 

It is instructive to compare (|4.12|) with an expression obtained by Baacke and Lavreshavili 
I22I ] for the same quantity (note that the collective coordinate term, —2 In (f^), is separated 
out in j^])' & l so using the MS renormalization scheme : 

Tbl = iln|r (P) (o )|+4iln7j 1) (oo) + ^< ) -^) 

i 00 r r i 21 1 

+ 2E( / + 1 ) 2 InfeHj-^H- ^(00) --(^(oo)) (4.15) 

The second term corresponds to the fourfold degenerate zero modes, and is computed in 
22l | by evaluating the slope at zero with an additional parameter added to the operator in 
order to displace the determinant from zero. This agrees numerically with our expression 
in (J3.13)) . [allowing for the collective coordinate term, —2 In (f^-)], but is more cumbersome 
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to evaluate. The subtractions h\ (op) and /ip^oo) in ()4.15|) are associated with first and 
second order Feynman diagrams 22j, and are given explicitly as the asymptotic values of 
h\ i (r) and hf\r), satisfying the differential equations 

1 



h\ 



(i)u 



(2)// 



■JZ+ll 



+ 2-— 



hU'-V(r) 



V{r)h\ 



(i) 








(4.16) 



with initial value boundary conditions: hi(0) = 0, and h'^O) = 0. The remaining terms in 
fl4.15j) are Born approximation terms 

1 



.4 



(i) 
fin 



drr V(r) 



A 



(2) 
fin 



128tH 



dq q 3 



V(q) 2 - vV + 4 In 



yV + 4 + q 
^fTl-q 



(4.17) 



where V(q) is the Fourier transform of the fluctuation potential V(r). 

Note that Baacke and Lavrelashvili's expression ()4.15|) also gives a finite answer for the 
renormalized effective action. This means that we have another, quite different, expression 
for the large I behavior of In (T(/)(oo)): 

1 



In (T®{oo)) ~ hPioo 



h 



(2), 



oo 



oo. 



+ 



(4.18) 



Comparing with ()4.14|) . we deduce that the difference between these large I behaviors should 



agree to order O y J • We have verified this numerically, and in fact we find that there 
is a nonzero difference at O 



/i, (1) (oo) + 



oo 



\$™drrV{r) |/ °° dr r 3 V(V + 2) 



O 



(4.19) 



(Z + l) (Z + l) 3 J V(^ + l) 5 . 

Similarly, comparing the terms outside the I sums in f!4.12|) and ()4.15|) . we see that these 
terms in (J4.12j) are local in the fluctuation potential V(r), while Ag^ in (|4.15|) is nonlocal in 
V{r). Nevertheless, the total quantities Tjjg and Ybl are equal (allowing for the collective 
coordinate term, — 21n(|^), that is separated out in j^)? and we have confirmed this 
equality numerically. This means that ()4.12|) and ()4.15|) are different ways of regulating the 



I sum, and since each is renormalized using MS, they should indeed be equal. In other 
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words, in (j4.12|) and (j4.15|) . the different subtractions inside the I sums are compensated 
for by different finite pieces outside the I sums, in such a way that the total agrees. This 
is very similar to behavior found in the computation of mass shifts in soliton theories (see 
Appendix B of 34]). It is also a highly nontrivial check on both expressions. But, although 
the expressions (J4.12|) and (J4.15j) are equal, in purely pragmatic computational terms the 
local expression (j4.12j) is considerably simpler. First, the zero mode part is more easily 
and more accurately computed using (j3.13|) . since it only requires the asymptotic properties 
of the (already computed) bounce solution $ c i(r). Second, one does not need to solve 
the differential equations ()4.1(ij) to find h\ (oo) and h\ ; (oo); instead, the large I behavior 
of lnT(£)(oo) is given by local expressions in V(r). Third, one does not need to compute 
(numerically) the Fourier transform of the fluctuation potential, nor integrate the result over 
all momentum, as is necessary to compute A-i in (j4.17|) . Instead, all subtraction terms in 
(I4.12J) are purely local in the fluctuation potential V(r). 

In a series of papers studying the thin-wall approximation, Konoplich and Rubin 

HQ 

have computed this same quantity with a related MS renormalization prescription, and 
obtained the following result in the a — > 1 limit [we have translated their notation to match 
oursl : 



9 

KR 



32 



i 2 * 



9V3. 



L — a) 

Their renormalization prescription is such that this leading term can be compared directly 
to ours. The subleading terms in ()4.20j) are scheme dependent, except for a universal log- 



arithmic term 



111 Q]. The agreement with our result (j4.12j) is striking. Because of the 
divergence as a — > 1, for plotting purposes it is useful to extract the overall factor of Tjj^p ■ 
In Figure |H1 we plot our result 1)4.12)1 . rescaled by a factor of (1 — a) 3 , compared with the 
Konoplich- Rubin thin- wall approximation answer ()4.20j) . also rescaled by a factor of (1 — a) 3 . 
The intercept at a = 1 is given by (j4.20j) to be 

— 1-— ^= ^ 0.16788888 , (4.21) 
32 I V ; 

which matches perfectly the a — > 1 limit of our result (|4.12|) . as can be clearly seen in Figure 

El 

The MS prescription provides simple renormalization structures, but we can also use a 
"physical" renormalization prescription. For example, we may impose the following normal- 
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FIG. 8: Plot of the MS renormalized effective action, I^jg, in ()4.11j) . rescaled by a factor of (1 — a) 3 , 
as a function of a in the thin- wall approximation limit as a — > 1. The solid dots denote our data 
points for a in steps of 0.01 from 0.9 to 0.99, and the solid curve is an interpolation through these 
points. The horizontal line is the leading term of the Konoplich- Rubin result (|4.20|) . also rescaled 
by a factor of (1 — a) 3 . Note the excellent agreement at a = 1. 



ization condition to the one- loop effective potential, 

d 2 X 
— {U cl + U eS ) U =(t>l = -(30^ - a 2 ), 



d 4 



(U d + U eS ) 



3A. 



In (I4.22J) . 0i = + 5(f)- is the solution of the one-loop modified field equation 



4i (Ud + U cS ) 



0. 



(4.22) 



(4.23) 



These normalization conditions can be satisfied when we change the renormalized parameters 
by finite amounts as A + 5X — >• A + 5X + Sf in X, and Aa 2 + 5(Xa 2 ) — > Aa 2 + 5(Xa 2 ) + 5f in (Xa 2 ). 



V. CONCLUSIONS 

In conclusion, we have presented a simple and practical technique for evaluating the 
prefactor determinant in the expression for the metastable decay rate in scalar field theories. 
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The technique for computing the determinant is based on the Gelfand-Yaglom formula 
for computing the determinant of a one-dimensional (here, radial) differential operator in 
terms of the asymptotic boundary value of an associated differential equation with initial 
value boundary conditions. This technique is extremely easy to implement for any given 
partial wave, but the naive sum over partial waves is divergent. Thus the direct application of 
the Gelfand-Yaglom formula to higher dimensions is not possible. However, this divergence 
can be regulated in various ways. Here we propose using the angular momentum cutoff 
regularization and renormalization scheme, which has been used previously to compute the 
explicit mass dependence of the fermion determinant in QCD for massive quarks in an 
instanton background 21]. In this approach the contribution of the low partial waves is 
computed numerically using the Gelfand-Yaglom formula, while the contribution of the high 
partial waves is computed analytically using radial WKB, which is a good approximation for 
large I. The merging of these two parts involves renormalization and we have illustrated our 
iechnique using an MS scheme in order to compare with the work of Baacke and Lavrelashvili 
^ , who also use the Gelfand-Yaglom technique, also with an MS prescription, but with a 
different regularization technique. We have shown that the two techniques agree, and have 
argued that the one presented here is somewhat easier to implement as it is purely local in 
the numerically determined fluctuation potential. The conversion to other renormalization 
shemes can be done using conventional field theory techniques, and corresponds to including 
finite polynomial terms in the effective potential. In an appendix we have derived a simple 
new formula for the / = 1 determinant with the zero modes removed, solely in terms of the 
asymptotic values of the bounce solution. 

The goal of this work has been to reduce the computation of field theoretic one loop 
fluctuation determinants to a straightforward numerical exercise. As a by-product, we have 
found in ()4.12|) a simple extension of the Gelfand-Yaglom result to higher dimensional radi- 
ally separable problems. There are many possible further applications of this technique, as 
there are many semiclassical problems where the classical solution, about which one is com- 
puting the quantum fluctuations, has radial symmetry. Closely related possible applications 
include: (i) metastable decay in theories with more than one field jj^j], where analytic and 
approximate approaches to the pref actor are quite difficult, but a direct numerical approach 
might be more useful; (ii) models in dimensions other than 4, which have been studied in and 



beyond the thin- wall approximation jfj, UjJ UJl LlZJ] ; (iii) an extension to finite temperature, 
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where the high temperature limit is essentially a dimensionally reduced 3d radial problem 
Q; but for intermediate temperatures the explicit summation over Matsubara modes 
is necessary. This technique for computing precisely the fluctuation contribution may also 
be useful for a set of fascinating questions concerning the validity of Lang er's homogeneous 
nucleation picture itself, as well as the semiclassical approximation 38[. More generally, 
our technique for extending the Gelfand-Yaglom formula to higher dimensions can also be 
applied to other symmetric semiclassical configurations such as vortices and monopoles (in- 
stantons were considered already in j2l|). 
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VI. APPENDIX A : ZERO MODE CONTRIBUTION 

In this Appendix we present a derivation of the expression (j3.13|) for the factor contributed 
by the / = 1 zero modes. We adapt a method of McKane and Tarlie 



fsee also 



28 



29 



40||). First, add a small quantity, k 2 , to the operator possessing the zero mode, and to the 
corresponding free operator (although this latter addition is not important in the end). Then 
for small k 2 the existence of the four zero modes for Ai [we suppress the (I = 1) subscript] 



implies that 



det(X + F) det^j 



det(M hee + k 2 ) detA^ free 
The result f!3.4|) means we can evaluate the LHS of ()6.1|) for arbitrary but small k 2 as 

_f (^mv , (6 . 2) 



det(A4 free + k 2 ) \r^oo \^ 2 ee {R) . 
where ip^ satisfies 

-< 2 -^ 2 + (^ + l + k 2 + V(r)^ k2 = , 

<M0)=0 ; ^(0) = 1. (6.3) 

Note that for / = 1 we can express the r — > boundary conditions (J3.6|) in this initial value 
form. The function ^ r 2 ee satisfies the same equation and boundary conditions as in ()6.3|) . but 
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with the potential V set to 0. So, one possible approach, as suggested 



in Q, 



to computing 



det' is to compute the derivative of lim^oo ( jjpssj^ j as k 2 — * 0. But a more direct and 
accurate method is as follows. 

Define the function ipo(r) to be the solution of (jfi.Hj) with k 2 = 0. Then by elementary 
integration by parts it follows that 



= J drr^ j-i (rV fca )' + + 1 + k 2 + V(r)^ 



[r 3 Oo?/V - ^o^jfea)] + ^ 2 / drr 3 ip ip k 2 . 



R 



(6.4) 



Applying the boundary conditions at r = we obtain 



MR) 
[MR) 



M(R) 



k 



2 dr r 3 -0o^ fc2 
(-J2Vo(fl)) 



(6.5) 



_ipk*(R). 

The important observation now is that V ; o( r ) is actually the normalizable zero mode of Ai, 
and decreases exponentially as e~ r at large r. On the other hand, for arbitrarily small 
but nonzero k 2 , the solution ipk 2 { r ) increases exponentially as e rv/1+fc2 at large r. Thus at 
large R and arbitrarily small but nonzero k 2 , the identity ()6.5|) implies that the leading k 2 
dependence at small k 2 is : 



k 







(6.6) 



Then ()6.2)1 leads to the following expression for the determinant with zero mode removed: 



det'A* 
det.M frcc 



J °° dr r 3 ipQ 



lim^oo [-2i? 3 ^(i?X ee ( J R)] 



(6.7) 



To simplify this general expression further, recall that the zero mode ipo can be expressed 
in terms of the classical bounce solution : 



M r ) 



*£(o) 



(6.8) 



This has two important implications. First, the numerator on the RHS of ()6.7j) can be 
expressed in terms of the classical bounce action : 

1 



dr r 3 ip 2 



2vr 2 (<&£(0))' 
1 

2vr 2 (^(O)) 1 



(6.9) 



23 



Since g^rbjjssg ^ s raised to the power —1/2 on the LHS of (|3.13|) . we see that the Sd^d] 
factors cancel. The second important implication of ()6.8|) is that the leading large R be- 
havior of ip' Q (R) is determined by the coefficient $oo in ()2.5|) . Since V'o ( r ) = ^(r)/?", the 
denominator in ()6.7)1 is 

$00 /A',(/-) 



{-2R 3 ^ (R)^ e (R)} ~ -2i? 



$ d (0) 



*£(o) 



8/ 2 (^) 

#^00 . (6.10) 



Thus we obtain the simple formula 



-1/2 



The final result ()3.13j) follows by noting that $"[(0) may be expressed in terms of $ c i(0) 
using the bounce differential equation ()2.1)1 and the boundary condition (|2.2)1 . 

We stress that the final expression ()3.13|) only requires knowing the two asymptotic 
constants, and $ c i(0), each of which is determined to great precision in the course of 
finding the classical bounce solution numerically, as described in Section IHII 
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